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ABSTRACT 


In this report. Part 2 of the final report, a perturbation model 
to the Marshall Space Flight Center (MSFC) Global Reference Atmosphere 
Model (GRAM) is developed for use in Aeroassist Orbital Transfer Vehicle 
(AOTV) trajectory analysis. The model reflects NASA Space Shuttle ex- 
perience over the first twelve (12) entry flights. The GRAM was selected 
over the Air Force 1978 Reference Model because of its more general for- 
mulation and wider use throughout the NASA. The add-on model, a simple 
scaling with altitude to reflect density structure encountered by the 
Shuttle Orbiter (shears, waves, "potholes-in-the-sky", etc) was selected 
principally to simplify implementation. Perturbations, by season, can 
be utilized to minimize the number of required simulations, however, 
exact Shuttle flight history can be exercised using the same model if 
desired. Such a perturbation model, though not meteorologically motivated, 
enables inclusion of High Resolution Accelerometer Package (HiRAP) re- 
sults in the thermosphere. Provision is made to incorporate differing 
perturbations during the AOTV entry and exit phases of the aero-assist 
maneuver to account for trajectory displacement (geographic) along the 
ground track. 



I . Introductory Background 


Development of a shear model for AOTV trajectory use based on 
Shuttle derived atmospheric data completes the requirements under 
the subject Contract. The model must include (at a minimum) the 
shear amplitude and frequency content which has already been observed 
in the STS accelerometry during the various entry flights. Yet, even 
though there is a large data base of Shuttle flights, there are no 
assurances that the "worse case" atmosphere has been encountered. 
Further, each flight exhibits (somewhat) its own unique density sig- 
nature. Thus, development of a simple, all encompassing, model is 
nontrivial, requiring, at least, the method of omphaloskepsis to begin. 
However, in support of AOTV activity, a model can be developed, and 
taking advantage of some seasonal similarities, the number of simula- 
tions can be minimized. This Section summarizes some considerations 
necessary in developing such a model. Surely there are others but it 
would seem that the major points to consider are the following: 


and, 


• utilization of STS specific atmospheres if desired, 

• selection of a baseline nominal (comprehensive) 
atmospheric model, 

• implementation, 

• minimization, 

• interpretation. 


Part 1 of this final report presented comparisons of Shuttle der- 
ived atmospheric parameters based on the first twelve STS flights. 
Comparisons were presented between the Shuttle derived density (and 
temperature) with those estimates provided by the National Weather 
Service remote sounding information as well as the MSFC GRAM and Air 
Force 1978 Reference Atmosphere models. The extent of density structure 
encountered during the three year STS flight history was quantified. 

It was stated that atmospheric perturbation models to emulate STS exper- 
ience resulted in significant AOTV trajectory departures, specifically 
the simulations reported as Reference 1 therein. Whether or not reason- 
able models have been (or can be) developed requires comparisons of the 
simulated dispersions with those actually encountered by Shuttle. 
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Specific Shuttle atmospheres can, of course, be developed as was 
done for STS-2, 4, 6, and 7. The method utilized was to develop a scale 
factor (using the 1962 Standard Atmosphere as reference) to operate on 
the modelled atmosphere to produce the equivalent density profile. Re- 
sults based on atmospheres generated in this manner were presented in 
Reference 1 of this report. Given the number of STS flights currently 
available, and considering the fact that this data base is rapidly in- 
creasing, continued evaluations of Shuttle specific atmospheres would be 
quite cumbersome, albeit expensive. Thus, where possible, it is advan- 
tageous to simplify the modelling requirements and still reflect the basic 
Shuttle environs encountered. This should not be construed as first 
priority but is certainly worthy of consideration. Other factors to 
consider, of course, are ease of implementation as well as interpreta- 
tion of the resultant perturbations modelled. For these reasons, though 
perhaps subject to argument, it is felt that the meteorologically motivated 
models of Robertson (Reference 2) might prove to be cumbersome. Finally, 
as alluded to in Part 1, there exists acceleration data which enable 
extension of the results above the threshold of the Inertial Measurement 
Unit (IMU) measurements, i.e., the HiRAP yg source can be utilized into 
the thermosphere. A model developed which can readily be expanded to 
higher altitudes to include the HiRAP results has some merit. 

Given the above considerations, a reasonable model can be developed 
for AOTV application. One needs to select a reference model and, to 
minimize the "degrees-of-freedom", reduce (reasonably) the STS atmos- 
pheric data base. Based on the results of Part 1, with the exception 
of the September flights and, to a certain degree on STS-9, the GRAM 
model appears to be the best reference data. This model has the added 
advantage that, unlike the AF'78 model which is only defined up to 90 km, 
it has a Jacchia-Roberts formulation for the uppermost altitudes. The 
GRAM also is more widely used throughout the NASA and, if required, has 
a spherical harmonic wind model available. With the GRAM as a baseline, 
and recalling that the STS flight derived atmospheres were (reasonably) 
similar by season, the limited perturbation model is next developed. 

In the following Section it is assumed that AOTV requirements , 
due to ground track considerations during the aero-assist portion of 
flight , will optionally want to include different perturbations during 
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the entry and exit phases of the maneuver. It would be remiss not to 
provide for this capability in view of the apparent "local" structure 
sensed by the STS accelerometers . Also , perturbing the nominal density 
from any source will not satisfy the gas laws. One might want to re- 
solve, say for temperature as an example , to preserve this relationship. 
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II . Proposed Perturbation Model 

Appendix A of this report shows Shuttle derived densities for each 
of the twelve (12) flights normalized to the GRAM values. Straight line 
segments are drawn through the curves to depict what is considered to 
be significant density structure.* As was done in Part 1, no flight 
updates to the predicted normal force coefficient were utilized and, as 
a consequence, each derived density curve could be shifted (more 
dense) by 3 to 5 percent to reflect the current Flight Assessment Deltas 
(FADS). The shaded region on these twelve figures is the GRAM uncer- 
tainty (±la) generated by the model. These are the MSFC data and not 
the computed statistics presented in Part 1, though reasonably compar- 
able. These data would be those utilized by analysts when performing 
Monte Carlo simulations. These data are superimposed for information 
only but, as a side activity, could be reviewed to characterize the 
adequacy of the GRAM statistical model in this altitude interval. 
Typically, each flight shows an approximate 2a departure from the GRAM 
in certain (restricted) intervals. These deviations are systematic, 
not random, occurrences. This might be considered adequate, however, 
comments would be appreciated. 

Figures 1 through 3 show composite plots of the results by seasons 
for spring, summer, and fall, respectively. Winter results can be seen 
(on the basis of a single flight) in Figure A-10 for STS-11. One can 
detect the density similarities from these seasonal composite charts, at 
least in the spring and summer months. The fall data presented would 
also show more similarity if STS-9 were disregarded. This flight is the 
outlier principally due to latitudinal effects. 

With caution, one could select the scaling results from a given 
flight for each season and, by providing for a datum shift, reasonably 
replicate the other flight results. This is, of course, not exact but 
worth considering to minimize the number of Shuttle perturbation models 
required, keeping in mind the (potential) factor of two to model separate 

^Specific STS flight results can be modelled by inputting a table of 
scale factors versus altitude whose break-point altitudes are as shown 
thereon. Certainly, much less than 50 pairs [hi, (Pcn/ p GRAM) J would 
be required. The resulting density would be the product of the scale 
factor (linearly interpolated with altitude) times the GRAM value and 
would thereby reflect the suggested shear structure. 
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entry and exit phases. Typically, in the spring months, STS-1 scaling 
(read directly from Figure A-l as an example) could be utilized. By 
shifting these results, for example by ± 20 kft, much of the signal 
shown in the remaining three flights would be represented. Similarly, 
in the summer months, the STS-4 results (as an example) could provide 
the benchmark perturbation model which could be shifted to emulate the 
others. One should keep in mind that the two September flights (STS-8 
and STS- 14) can be shifted by approximately 20 percent at the top due 
to the fact that the GRAM is too dense in this month. In fact, the 
results can be essentially rotated about h~230 kft, wherein the excess 
model density first becomes apparent (see Part 1), peaking to ~20 percent 
at h~320 kft. Such an adjustment would make these two flights fall more 
in line with the other two summer flights. It is noted that the sugges- 
ted shear structure in the narrow vicinity of 230 kft to 250 kft for the 
summer months is very repeatable. Shifting the STS-4 data (or the re- 
sults from any of these flights) by ± 20 kft must be reviewed. If that 
specific structure, when displaced, creates a significant problem for 
A0TV analysts, the results must be carefully reviewed. Lastly, either 
the STS- 2 or STS- 5 values can be used for the fall scaling with, again, 
a ± 20 kft datum shift provided. This, of course, disregards the STS-9 
flight results. 

It is recommended that a representative flight be selected (each 
analyst might want to make his or her choice) and simulations be con- 
ducted in which the scaling versus altitude is shifted as suggested 
(±20 kft). Those simulations which show potential problems must be re- 
viewed versus actual historical STS flight data. The charts herein in 
Appendix A can serve that purpose. Future flights, if more radical 
density structure is encountered, will need be factored into the analysis. 
As the AOTV design solidifies, and dependence on the various density 
structure (magnitude, altitude of occurrence, etc.) is better established, 
a further review of actual STS atmospheric phenomena might be warranted. 

To implement this model a table(s) of scale factors versus altitude, 
a linear interpolation scheme to operate on the table(s), as well as a bias 
value (s) on altitude (to shift the table (s)) would be required inputs, 
assuming of course the GRAM is already available in the user's software. 
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If not, STS atmospheres scaled to the '76 Standard (see Part 1) could 
just as readily be utilized. One could implement a single table per 
season and adjust the altitude datum bias through some mechanism based 
on AOTV altitude rate to model differing perturbations during the entry 
and exit phases. In any event, trivial software modifications would be 
required to implement the following schematic: 



Readers are reminded that the charts of Appendix A are ratioed to 
the comprehensive GRAM values, to include latitudinal and seasonal 
effects. If the 13 th month "average" model is utilized the scale fac- 
tors shorn, should still be adequate for analysis purposes during the 
preliminary design stages , e.g to evaluate and quantify the effects 
of "representative" STS density structure on AOTV trajectories. Also, 
in addition to indicating structure, the scale factors reflect syste- 
matic, local, deviations from the GRAM estimates. These curves (or 
perhaps a smooth version of same) could therefore be utilized as a 
basic shift to the GRAM data in conjunction with the GRAM statistics 
in a Monte Carlo analysis . Thus the typical "net" zero mean error re- 
sulting from the Monte Carlo runs would actually end up describing, 
apart from the random signal, the baseline offset to the GRAM evidenced 
in the STS data. 
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III. Extension to Include Perturbations in Thermosphere 

Figure 4 shows STS-6 and STS-7 HiRAP results (see Reference 8 of 

Part 1) obtained from R. C. Blanchard of the Aerothermodynamics Branch, 

Space Systems Division of LaRC. With yg resolution, derived densities 

well up into the thermosphere are obtainable. The data, shown as a 

straight line segment of Blanchard's results, are contiguous with the 

IMU derived values as one would expect. HiRAP saturation occurs at 

8 mg's wherein the IMU signal, though noisy, becomes meaningful. Shown 

on each sub-figure of Figure 4 are actually two curves; 1) Blanchard's 

results scaled to the '62 Standard, and 2) the GRAM scaled to the same 
2 

Standard. The HiRAP data show large waves in the thermosphere, similar 
in nature to that evidenced in the GRAM (perhaps substantiating in part 
the Jacchia-Roberts formulation) but much larger in amplitude and, in 
the case of STS-7, shifted by some 50 kft in altitude. Figure 5 shows 
the same two flights normalized to the GRAM. It is recognized that this 
is an extremely high altitude, low density, region but it is not incon- 
ceivable that such density departures would have some effect on AOTV 
trajectories, particularly during the exit phase. Thus, with virtually 
no a priori experience, it is recommended that consideration be given to 
implementing the data of Figure 5, appended to the upper end of some 
scaling table, as part of the AOTV trajectory analyses and design. 


2 

Design values were utilized in the GRAM for actual and mean solar 
flux as well as for geomagnetic index. 
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IV. Conclusions 


An atmospheric perturbation model for AOTV analysts has been 
developed to enable trajectory simulations using density profiles 
commensurate with those sensed during the STS Orbiter entry flights. 
Attempts have been made to minimize the requirements by evoking sea- 
sonal similarities observed in the Shuttle data. The proposed model, 
which should be readily implemented and the resultant perturbations 
easily interpreted, can be utilized during the entry and exit phases 
of the aero-assist maneuver portion of flight. Actual conclusions 
as to the adequacy of the model require user feedback as well as con- 
tinued investigations of atmospheric structure encountered on future 
STS flights. For the most part, a sizeable data base of Shuttle results 
is available and, it is felt, the essence of the encountered atmos- 
pheric structure is adequately modelled for AOTV design studies. 
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APPENDIX A 


Comparison of Shuttle Derived Densities 
with MSFC GRAM Data for 
First Twelve Shuttle Entries 

Notes : 

• No FADS utilized in derivation 

• Shaded region shown is ±la band as 
defined by MSFC for their model. 
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Figure a-5. STS-5 (November) GRAM density scaling 
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Figure A-io. STS— 1 1 (February) GRAM density scaling 
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